Anomalous scaling at the quantum critical point 
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We show that Hertz cj> 4 theory of quantum criticality is incomplete as it misses anomalous non- 
local contributions to the interaction vertices. For antiferromagnetic quantum transitions, we found 
that the theory is renormalizable only if the dynamical exponent z — 2. The upper critical dimension 
is still d — 4 — z — 2, however the number of marginal vertices at d = 2 is infinite. As a result, the 
theory has a finite anomalous exponent already at the upper critical dimension. We show that for 
d < 2 the Gaussian fixed point splits into two non-Gaussian fixed points. For both fixed points, the 
dynamical exponent remains z = 2. 



Quantum phase transitions (QPT) at zero tempera- 
ture are currently subject of intensive experimental and 
theoretical study. These transitions occur in a number of 
itinerant fermionic systems under the change of pressure, 
doping, magnetic field, or some other external parame- 
ter. QPT are very different from finite temperature phase 
transitions as the dynamics of the order parameter field 
can be neglected at finite T, but cannot be neglected at 
T = 0. The conventional phase diagram in the (x, T) 
plane, where x is the external parameter, has three dis- 
tinctive areas (see Fig. P). The ordered phase (for most 
cases, antiferromagnetic) is to the left of the quantum 
critical point (QCP), the disordered, paramagnetic Fermi 
liquid phase is to the right, and right above QCP there 
is a quantum critical regime that we will study. The sys- 
tem properties in this regime are governed by quantum 
dynamics of slow fluctuations of the order parameter . 

In his original approach to quantum criticality for itin- 
erant fermions, Hertz Q considered coupling between 
fermions and low-energy bosonic field which condenses 
at QCP. He integrated out fast fermions and obtained an 
effective theory for the slow bosonic degrees of freedom 
which, he argued, is Landau-Ginsburg- Wilson (LGW) </> 4 
theory with the upper critical dimension d cr = 4 — z, 
where z is the dynamical exponent. At d = d cr , the 
<j) A vertex is marginal, above d cr it is irrelevant, and for 
d < d cr it is relevant. Higher order if) 6 , etc vertices are 
all irrelevant near d cr . This theory was latter extended 
by A. Millis Q and others Q to explain the finite tem- 
perature properties of metals in the vicinity of QCP. 

In recent years, the applicability of Hertz-Millis theory 
to quantum phase transitions in heavy fermion metals 
has been questioned. The theory seems to work in some 
systems and do not work in others (see, e.g. It has 
been suggested that this inconsistency may be due to the 
fact that QPT in some heavy fermion materials require 
two-band description and may be accompanied by the 
discontinuous change of the area of the Fermi surface . 
This phenomenon is not included in Hertz theory which 
describes one-band itinerant models. 



In this paper we consider QPT for which one-band de- 
scription is valid, and the Fermi surface changes continu- 
ously through the transition. We show that even for these 
transitions, Hertz's theory is incomplete. Specifically, we 
argue that upon integrating out fermions, the effective 
bosonic theory becomes non-local. We show that only 
the theory with z — 2 is renormalizable. For d > d cr = 2 
all interaction terms are irrelevant, and the fixed point is 
just Gaussian. However, for z = 2 and d = 2, there is an 
infinite number of marginal terms in the effective action. 
For d < 2, the Gaussian fixed point splits into two new 
fixed points, one stable and one unstable. The scaling 
dimensions for different components of the momentum 
are different for these new fixed points, and for both of 
them, the dynamical exponent z measured in units of the 
most relevant component of the momentum equals 2. 

The existence of an infinite number of marginal terms 
at d = 2 is in variance with the LGW theory in which 
only <^> 4 vertex is marginal. A single marginal vertex gives 
rise to only logarithmic corrections to a Gaussian theory 
and does not change critical exponents. This allows one 
to expand critical exponents in e = 2 — d. An infinite 
number of marginal vertices, however, gives rise to much 
stronger, power-law corrections in which case e expansion 
does not work, and anomalous exponents emerge already 
at the upper critical dimension. This, in particular, ex- 
plains why explicit computations of the dynamical spin 
susceptibility |j| at the antiferromagnetic QCP in d = 2 
yield x(^,q) K ((Q — q) 2 + |^m|)~ 1+7 with the anoma- 
lous exponent 7 ~ 0.25. 

Our point of departure is the same as in Hertz theory 
- we consider fermions coupled to a bosonic field cf> whose 
dynamics is governed by the dynamical exponent z. The 
Lagrangian density of this model has a form 

L = c w , k G ( 7 1 (w,k)c a) , k + ~ XoH^q^n.q + 

gC Wik C w _SXk-q0n,q (1) 

where c^,k, and c Wj k are Grassmann variables, g is the 
coupling constant, and the bare fermionic and bosonic 
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FIG. 1: The schematic phase diagram for Quantum Phase 
Transition. 





FIG. 2: The diagram for the <j> term in the effective bosonic 
action (left) and the schematic diagram for <f> 2n (right). Dif- 
ferent colors show Green's functions near different hot spots. 



propagators are 

Gq 1 (w,k) = uj - e k 
Xo 1 (r»,q)=r 2 + (q-Q) 2 + (^) 
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Here ek is the quasiparticle excitation energy, £ is the 
correlation length which diverges at QCP, and Q is the 
ordering vector for cf>. An input for the theory is the 
assumption that there exist special hot spots (or lines) 
at the Fermi surface separated by Q. In the absence of 
nesting, the Fermi velocities at and k/ ls + Q are not 
antiparallel. 

The action of Eq. is quadratic in fermions, and 
Hertz proceeded by integrating fermions out. This renor- 
malizes Xo(fi,q) an d also introduces the interaction be- 
tween bosonic fields: <fi 4 ,(f> 6 terms, etc. Due to particle- 
hole symmetry, only terms with even number of bosonic 
field are generated. The resulting action is 8] 



s H = i fdoAx" l («,q) 

f; f (dnd 2 ?) 2 "- 1 ^ 
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The renormalization of Q) an d the vertices bi n 
are given by the diagrams shown in Fig. [2] Hertz as- 
sumption that the vertices are local implies that they all 
can be evaluated at bosonic frequencies equal to zero, 
and momenta equal to Q. Consider & 4 as an example. 
The corresponding diagram contains two pairs of Green's 
functions with momenta near k; ls and k; ls + Q (see Fig. 
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One can easily make sure that once we linearize the 
dispersion near hot spots, the integral over momentum 
would vanish because of double poles 3]. This implies 
that the integral in (|SJ) comes from electrons with high 



energies, of the order of bandwidth W, where the spec- 
trum cannot be linearized. Accordingly, the value of 64 
scales inversely with W and vanishes in the continuum 
limit W — > 00. The same consideration holds for all other 
&2n- For finite W, all &2n ar e just some constants. It is 
then straightforward to calculate the engineering dimen- 
sionalities of the couplings. Using the fact that [SI] = z 
(see Eq. ©), we immediately find that in d = 2 the 
dimensionality of the field is [4^ ] = — 4 — z, and the 
dimensionalities of the vertices are [&2n] = 2 — (n — l)z. 
For z = 2, this implies that 64 is marginal, while all the 
other vertices are irrelevant. This is the known result of 
the Hertz theory HQ. 

We now demonstrate that the assumption of the lo- 
cality of the interaction is in fact incorrect. Indeed, we 
found above that the integral in J5j is zero for a linearized 
dispersion, because of double poles. However, this is true 
only when the bosonic frequencies Sli are zero. If we con- 
sider instead the limit Sli — ► 0, we find that there exists 
a tiny range of \u>\ < \Sli\ where the double poles split 
into pairs of closely located poles in different half planes. 
The momentum integration then results in small denom- 
inators. This gives rise to universal, anomalous contribu- 
tions to &2n which may be quite large if the small denom- 
inator overshoots the smallness of the frequency range. 

Linearizing fermionic dispersion near k/j S and k/ ls + Q 
and carrying out integration over momentum and fre- 
quency in Eq. (J^J for nonzero bosonic frequencies Sli and 
qi ^ Q, we find that 
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where q = q — Q. The exact expression has a complex de- 
pendence on all external momenta and frequencies 0, 
but since we are only interested in the engineering di- 
mensions, the estimate in (0 is sufficient. We see that 
64 strongly depends on the ratio Sl/vFq and actually be- 
comes large in the limit ri — > 0, vpq — 0. We emphasize 
that 64 in JSJ is universal in the sense that it does not 
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depend on the details of the fermionic dispersion at en- 
ergies of order W and survives in the limit W — * oo. 
Restricting with only universal piece in 64, we find that 
d> i term in the effective bosonic action becomes 



<? 4 / (dvnr 
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(7) 



where 34 ~ g i /v F . Performing the same calculation for 
vertices b 2n with n > 2, we obtain that 
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Accordingly, 
form 



6 2 ™ term in the effective action takes the 
101 
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, q ) 2 " (9) 



(0 - vpq + iS) 2 ^- 1 *) 



where 52™ oc g 2n /v F . 

We can now re-evaluate the scaling dimensionality of 
the vertices. Performing the same estimates as for the 
Hertz theory, we obtain 
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There are two consequences of this result. First, Eq. 110(1 
holds for all n, down to n = 1. For n = 1, the universal 
<?2 vertex generated by fermions is 



g 2 / (d 2 gdn)\n\ (^y 
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This is nothing but the Landau damping term. Adding 
the g 2 vertex to the Gaussian part of the bosonic action, 
we find that z = 2 is special in that the form of Xoi^i l) is 
reproduced, i.e., the bosonic dynamics is self-generated. 
Furthermore, z — 2 dynamics will obviously dominate 
even if in the bare Xo(^j q) z < 2. Second, we see from 
(|10(l that for z = 2 all vertices are marginal in d = 2. 
This is very different from the Hertz theory where only 
<j) A vertex was marginal. There, a single marginal vertex 
lead to only logarithmic corrections to xo(^, l) and did 
not modify critical exponents which remain mean field in 
d — 2. However, when the number of marginal vertices is 
infinite, this argument does not work as each vertex now 
gives rise to logarithmic corrections to the susceptibility. 
The infinite series of logarithms coming from all b 2n sum 
up into a power-law correction to the susceptibility, such 
that at £ = 00, 
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where 7 > 0. This implies that the system develops an 
anomalous exponent 7 already at the upper critical di- 
rl ' - -• 0.25_was earlier obtained 



mension. The exponent 7 — 

in perturbative 1/N calculations yj iui jzjq. _lih 
present consideration provides understanding of its ori 
gin. 



for Eq. ©. The 




FIG. 3: Hot manifolds in 3D. 

We next consider d ^ 2. For Hertz theory, simple 
power counting gives [4> 2 ] = —2 — d — z, and [&2n] = 
(2 — d — z)n + d + z. This implies that 



[64] = 4-d-z 



(13) 



For z — 2, 4 term is irrelevant for d > 2, marginal 
for d = 2, and relevant for d < 2. All other vertices 
are irrelevant near d — 2. However, as for d — 2, this 
power counting is invalidated because there exist univer- 
sal, non-analytic contributions to &2n which overshoot 
regular pieces. The calculation of the anomalous pieces 
in bi n for d ^ 2 is a bit more involved compared to d = 2 
case as in arbitrary dimension, there exist hot manifolds 
on the Fermi surface instead of hot spots (see Fig. |3J. 
Still, in a single scattering event, an electron is scattered 
from one hot manifold to another. We therefore can again 
split the Green's functions in the integrals for the vertices 
into the two groups belonging to one of two manifolds, 
and in each manifold evaluate the anomalous contribu- 
tion coming from momentum integration transverse to 
the Fermi surface. The integration over the remaining 
components of momentum gives the volume of the man- 
ifold. Performing this integration, and also integrating 
over frequency, we obtain 
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term in the action then becomes 



!/■,„ I (dV^) 2 "- 1 (fe q ) 2 " 

Simple power counting now yields 

[gin] =d-2-{d + z-4)n 



(14) 



(15) 



(16) 



We see that the Gaussian part [52] = 2 — z independent 
of d, i.e., z — 2 dynamics is self generated in any dimen- 
sion. We also see that for z — 2, g 2n — — [n — l)(d — 2). 



4 



This implies that for d > 2 all vertices with n > 1 are 
irrelevant, for d — 2 they are all marginal, and for d < 2 
they are all relevant. We see therefore that for d < 2, the 
number of relevant vertices is infinite. This is is another 
discrepancy with the Hertz theory. 

The case d < 2 requires further study. Previous analy- 
sis of the scaling dimensions of the vertices was pertained 
to a Gaussian fixed point. Below d = 2, this fixed point is 
no longer stable, and we need to find the new stable fixed 
point. Our strategy is the following: One can straight- 
forwardly verify that for d < 2, the integration over the 
hot manifold, leading to Eq. Ijl5|l . affects differently the 
components of the bosonic momenta parallel to the or- 
dering momentum Q and perpendicular to Q [q^ and 
qj_, respectively). We therefore introduce an extra scal- 
ing dimensionality r\ for one of them. Using r\ as input, 
we obtain scaling dimensionalities of the coupling con- 
stants g2n- We then use the fact that at the fixed point 
all vertices must be marginal, and obtain the values of z 
and 77 at the new fixed points. 

Let's measure the scaling dimensionalities in units of 
q«. Then, by definition, [qu] = 1, while for qj_ we intro- 
duce 
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The Gaussian part of the action is / (d d 1 q±dq\\dtl)(q 2 + 
q^)0 2 . For r]± > 1, we can neglect qj_ in comparison to 
<7H , while for t]±^ < 1 we can neglect g|| in comparison to 
q^. Using this, we find after a simple algebra that the 
dimensionality of the field is 



[(j) 2 } = -3 - (d - 1)77 -z, for 77 > 1 
[4> 2 ] = -1 - (d + 1)77 - z, for 77 < 1 



(18) 



Using Ijl8|) . we then find by straightforward calculations 
that the dimensionalities of the coupling constants gin 
are 

[<?2„]= n[3-(d-l)r ] -z] + [(d-l)r ) -l] 
for 77 > 1 

[S2n] - n[-l -(d- 5)77 -z] + [(d - 3)77 + 1] 

for 77 < 1 (19) 

(care has to be taken in regularizing the infrared diver- 
gences in the momentum integrals for 77 < 1). At the 
fixed point [<?2n] = for all n. Solving for [<?2n] = 0, we 
obtain 



v = i/(d-i) 

z = 2 

V = V(3 - d) 

z = 2/(3 -d) = 2?y 



for 77 > 1 
for 77 < 1 



(20) 
(21) 



exist two fixed points. They progressively deviate from 
each other as 2 — d increases. Obviously, one of these fixed 
points is stable and the other is unstable. To determine 
which one is stable, we use the result of perturbative 1/N 
analysis of logarithmic corrections to VFq\\ and VFq± for 
d = 2 0- According to Q, the ratio q±/q» oc 1/logf; 
renormalizes to zero at £ = 00 (this makes the Fermi 
surface nested at the hot spots at QCP). By continuity, 
this implies that the stable fixed point is the one with 
77 > 1. We see from (1211) that at a stable non-Gaussian 
fixed point the dynamical exponent remains z = 2 

In summary we have demonstrated that Hertz theory 
of quantum criticality is incomplete. We have shown that 
all vertices in the effective bosonic theory possess univer- 
sal, anomalous pieces which overshoot constant terms of 
the Hertz theory. For quadratic part of the action, this 
universal contribution is the Landau damping term. We 
found that z = 2 is special in that the theory is renormal- 
izable for any d. The upper critical dimension for z = 2 
is d cr — 2, as in Hertz theory, however, the number of 
marginal vertices at d — 2 is infinite. This gives rise 
to the appearance of the anomalous exponent already at 
d = 2. For d < 2, the number of relevant vertices is infi- 
nite. We found that the Gaussian fixed point is unstable, 
and there exist two non-Gaussian fixed points with dif- 
ferent scaling dimensions for different components of the 
momentum. At each of these new fixed points, there is 
an infinite number of marginal vertices. We found that 
one of these points is stable and one unstable. Still, for 
both fixed points, the dynamical exponent, measured in 
units of the largest component of the momentum, re- 
mains z = 2. 
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We see from (|20l I21f) that each inequality is indeed 
satisfied only if d < 2. This implies that for d < 2 there 
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